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Abstract 

Methods from scattering theory are introduced to analyze random Schrodinger 
operators in one dimension by applying a volume cutoff to the potential. The key 
ingredient is the Lifshitz-Krein spectral shift function, which is related to the scattering 
phase by the theorem of Birman and Krein. The spectral shift density is defined as the 
"thermodynamic limit" of the spectral shift function per unit length of the interaction 
region. This density is shown to be equal to the difference of the densities of states for 
the free and the interacting Hamiltonians. Based on this construction, we give a new 
proof of the Thouless formula. We provide a prescription how to obtain the Lyapunov 
exponent from the scattering matrix, which suggest a way how to extend this notion to 
the higher dimensional case. This prescription also allows a characterization of those 
energies which have vanishing Lyapunov exponent. 
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1 Introduction 

In this paper we consider random Schrodinger operators H(uj) in L^(]R) of the form 

i/H = /Jo + $^«,M/(--j), ^0 = -^, (1.1) 

where {aj{uj)}j^i is a sequence of i.i.d. (independent, identically distributed) variables 
having a common density Lp (i.e. P{aj G rfy} = Lp{y)dy), which is continuous and has 
support in the finite interval In what follows we always suppose that the single- 

site potential / is piecewise continuous with supp/ C [—1/2, 1/2]. Moreover, we require 
that / > 0. The spectral properties of ( |1.1| ) were studied in detail in [j3^, |16|, |31|]. The 
results are most complete for the case when / is the point interaction (see [Q]). 

The main idea of our approach is to approximate the operator H{uj) ( |1.1D by means of 
the sequence 

j=n 

H^''\uj) = Ho+J2 «.M/(- - j) (1-2) 

j=~n 

with unchanged Hq, which converges to H(uj) in the strong resolvent sense. This ap- 
proximation gives the opportunity to use scattering theory in order to study the spectral 
properties of the limiting operator (|1.1|). In fact, we show how to recover the spectral char- 
acteristics of H{uj) from the limiting behaviour of the spectral characteristics of H^"'^(uj) 
in the "large support" limit n ^ oo. 

One of the important ingredients of our approach is the Lifshitz-Krein spectral shift 
function (see [Q] for a review). Recently it has received renewed interest due to its appli- 
cations to different problems in the theory of Schrodinger operators [|2^, ^ ^ ^ |m ^ 
^ [T^, |2T|, ^ 0. In the context of our approach the spectral shift function 
naturally replaces the eigenvalue counting function, which is usually used to construct the 
density of states for the operator ([lT|). The celebrated Birman-Krein theorem [Q] relates 
the spectral shift function to scattering theory. In fact, up to a factor — tt^ it may be identi- 
fied with the scattering phase for the pair (if (cj), /fo), i-e. ^^''\E;uj) = -'k'^6^''\E;uj) 
when E > 0, 

d^^'HE-uj) = -logdet5(")(E;cu) = -logdet | ^ll^^^ 
^ ' ^ 2z ^ ' ^ 2« \^ L^^\e) 

Here |T(")(E)|2 and \R^''\E)\^ = |L(")(E)|2 have the meaning of transmission and re- 
flection coefficients, respectively, such that \T^''\E)\'^ + |i?(")(E)|2 = 1. For < 
^(")(E;u;) equals minus the counting function of H^"'(uj). 

These two properties of the spectral shift function, namely its relation to scattering 
theory and its replacement of the counting function in the presence of an absolutely contin- 
uous spectrum convinced the authors already some time ago that the spectral shift function 
could be applied to the theory of random Schrodinger operators. In fact, our previous arti- 
cles [ p4] , \[% ^ were in part preparatory investigations aiming at such an application. In 
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[|TPp we proved convexity and subadditivity properties of the spectral shift function with 
respect to the potential and the coupling constant, respectively. Such properties often show 
up when considering thermodynamic limits in statistical mechanics. In [0, R3p we studied 



cluster properties when the potential is a sum of two terms and the center of one is moved 
to infinity. Again such properties play an important role in statistical mechanics as well as 
in quantum field theory. Possible applications of the theory of the spectral shift function 
to random Schrodinger operators have also been envisaged by Simon [|58|]. Some other 
applications of the scattering theory in one dimension to the study of spectral properties of 



Schrodinger operators with periodic or random potentials can be found in [g9|, |5j] and [ plQ 
respectively. 

Below we prove some new inequalities for the spectral shift function, which reflect its 
"additivity" properties with respect to the potential being the sum of two terms with disjoint 
supports. These inequalities are closely related to the Aktosun factorization formula ( |2.25D 
[|T]] (see also [|5^, ^ |53|]). Combined with the superadditive ergodic theorem they will allow 
us to prove the almost sure existence of the limit 

m = lim VxT^' ^1-3) 

n^oo 2n + 1 

which we call the spectral shift density. We prove the equality ^{E) = Nq{E) — N{E), 
where N(E) and Nq(E) = 7r^^[max(0, -E)]^/^ are the integrated density of states of the 
Hamiltonians H(uj) and Hq, respectively. Also we reconsider the Aktosun factorization 
formula and show that it is a direct consequence of the propagator property of the funda- 
mental (or transfer) matrix of the Schrodinger equation. 

Another very important quantity associated with the Hamiltonian ( pTTj ) is the Lyapunov 
exponent 'j{E). In particular, according to the Ishii-Pastur-Kotani theorem [ |37| ] the set 
{E : 7(-E) = 0} is the essential support of the absolute continuous part of the spectral 
measure for H{uj) (|LT1). 

We wiU establish that for > the function -j{E) + m{N{E) - Nq{E)) can be 
interpreted as the density of the logarithm of the transmission amplitude T(E), 

\im^-^^^^ = -J{E)-^nm■ (1-4) 

n^oo 2n + 1 

The connection between the Lyapunov exponent and the transmission coefficient |Ti,"'* (E) \ 
was recognized long ago though a rigorous analysis was still absent. It is well 

known (see e.g. []r3[]) that the function w{E) = — 7(-E) + inN{E) can be analytically 
continued in the complex half -plane ImE > as a Nevanlinna function (i.e. an analytic 
function which maps the upper complex half plane into itself). We will recover this prop- 
erty of w{E) directly from the analytic properties of the transmission amplitudes T^"'\e). 
Moreover, as a direct consequence of these properties we give a new proof of the Thouless 
formula 

7(E) - jo{E) = - [\og\E- E'mE'), (1.5) 
Jr 

where 7o(-E) = [max(0, — i?)]^/^ is the Lyapunov exponent for Hq. 
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Also we prove a new representation for the Lyapunov exponent for positive energies, 



lim 

n-^oo 2n + 1 



where 



log 



n 



(1.6) 



Lc{E) 
TciE) 



RcjE) 
T^(E) 

tJey 



(1.7) 



and Ta{E), Ra{E), La{E) are elements of the S-matrix at energy E for the pair of Hamilto- 
nians (Hq + af, Hq). This representation allows us to apply the theory of random matrices 
to prove that 7(-E') > for a.e. > almost surely, which in turn by Ishii-Pastur- 
Kotani theorem implies that the spectrum of H has no absolute continuous component in 
(0, oo). We give also an explicit description of the set of special energies where 'y{E) = 0. 
To our best knowledge this set was known explicitly only for the two particular cases 
when the single-site potential / is a 5-potential or the characteristic function of the interval 
[-1/2,1/2] (see@). 

We express the density of states N{E) for positive energies in terms of the product of 
matrices 



N{E) = T- lim 7^— ^arg e±, TT A«^.H(^)e± 



(1.8) 



with Cj 



(1,0)^ 



(0, 1)^ and (., ■) being the inner product in C^. This representa- 



tion is similar to the definition of the density of states through the rotation number of the 
fundamental solution of the Schrodinger equation pE|]. 

Formulae (11.3D, (|1.4t), (11.6|), and ( ]1.8D also provide very simple and efficient numeri- 
cal algorithms to compute the density of states and the Lyapunov exponent of disordered 
systems in one dimension. In this context we also remark that the representations ( [L3| ) 
and ( |1.8[ ) for finite n and E > give smooth approximations to the spectral shift density 
and density of states, respectively. This contrasts with the usual procedure (see e.g. p8|]), 
where the density of states is approximated by step-like functions. As an example in Sec- 
tion 4 we have calculated the spectral shift density for the deterministic Kronig-Penney 
model. 

Our approach is also applicable to deterministic Hamiltonians of the form H = Hq + 
V with a potential V, which is supposed to be uniformly in ^('^^(M) and for technical 
reasons also bounded, but without any assumption on its decay at infinity. Let {yj}j^z 
be a sequence of real numbers such that yj ±oo as j ±oo and yj < yj+i for all 
j G Z. In this case we approximate H through H^"-^ = Hq + \/^{") with x*-"^ being the 
characteristic function of the interval [y_„, Again we show that the relations analogous 
to (O), (|13), dO]), (III) hold. 

An extension of ideas developed in the present paper to the case of higher dimensions 
will be given in [ p6| ] . 

Acknowledgements: We are indebted to J.M. Combes and V. Enss for valuable re- 
marks. 
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2 Auxiliary Results 

We start with a short discussion of some important properties of the spectral shift function 
^{E) = ^{E; H, Ho) for a pair of self-adjoint Hamiltonians H = Hq + V and Hq = 
—cP/dx^ on L^(]R) with domain of definition being the Sobolev space ly^'^(M) (see e.g. 

for the definition). Here V denotes the multiplication operator by the real valued 
function V{x), which is supposed to satisfy 

j {l + \x\'^)\V{x)\dx <oo. (2.1) 

The spectral shift function is defined by the trace formula 

tr (0(if) - 0(ifo)) = / <P\E)i{E)dE, (2.2) 

JR 



which is valid for a wide class of functions (see [^). For instance, elements in 
are in this class. Relation ( |2.2D defines i{E) only up to an additive constant. We fix this 
ambiguity by the condition that i{E) = for all E below the spectrum o-(H) of H. 

With this normalization condition the function ^ (E) for < equals minus the num- 
ber of the eigenvalues of H = Hq + V less than E, 

aE-0) = -nE-o{V). (2.3) 

For E > the relation between ^{E) and the scattering matrix S{E) is given by the 
celebrated Birman-Krein theorem [|^], 

log det S{E) = -2m^{E). (2.4) 

Here the branch of the logarithm is fixed by the condition ^(E) ^ as £" — > cxo (see 
Lemma 2.1 below). Since S{E) is continuous for all > [|T8|], so is ^{E). 

The relations ( [2. 3D and ( [2.4D define ^(E) everywhere, i.e. also for < 0, except at 



the finite number of points of discontinuity at the eigenvalues of H and possibly at = 0. 
We can redefine ^{E) by requiring that ^{E + 0) = ^{E) for all E G M such that ^{E) 
becomes right semicontinuous. 

We recall that the spectral shift function is monotone with respect to the perturbation, 
i.e. if Hi and H2 are self-adjoint operators such that Hi < H2 in the sense of quadratic 
forms, then ^{E; Hi, Hq) < ^{E; H2, Hq) (see e.g. 

The value of ^(0) depends on the spectral properties of the point E = 0. We call 
the point E = regular iff (/ + V^/^i^ol-z)!^!^/^)"^ exists and is bounded at 2; = 0. 
In the opposite case E = is called an exceptional point. Here we used the notation 
\V\^/\x) = |V(x)|i/2 and V^/^{x) = signV(x)|V(x)|i/2 ^nd Ro{z) = {Hq - z)'^ is the 
resolvent of the free Hamiltonian Hq. Other characterizations of the exceptional case can 
be found in [0, |8|, Ip. The Levinson theorem for Hamiltonians on a line 0, P states that 

m=a+0) = -no{V) + ^ (2.5) 
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if E = is a regular point for H and 

m = e(+0) = -no{V) (2.6) 



if = is an exceptional point. An extension of ( |2.6| ) for potentials with slower decrease 
can be found in [p6|]. 

The scattering matrix S{E) for the pair of Hamiltonians {H, Hq) at fixed energy E > 
is a 2 by 2 unitary matrix (see [ JTSl , p3|]) 

"^^^^ - ^ L{E) T{E) ) ■ (2-^^ 

Below we will use the fact that due to unitarity the S-matrix can be parameterized by the 
absolute value of the transmission amplitude < \T(E)\ < 1 and two real valued phases 
5{E) ande{E): 



V i^/l-\T{E)\''e^'^^)-^'^^) \T{E)\e^'^^) ) ' ^ ^ ^ 

Here is the scattering phase such that 5{E) = —n^(E). For reflection symmetric 
potentials exp{2i6{E)} = 1. 

Below we will need the following auxiliary results. 

Lemma 2.1 Let V{x) satisfy the condition f |2.ip . Then there is a constant Cy > such 
that 

\Cm<Cv (2.9) 
for all E . Moreover, there is a constant C > independent of V and E such that 



\^{E)\<C{^ / \Vix)\dx+ 



1 " 



\V{x)\dx 



(2.10) 



for all E>0. 



Proof. First we prove the inequality ( [2.101 ). By monotonicity of the spectral shift func- 
tion we have 

^E; Ho -\V\,Ho)< ^E) < ^E; Ho + \V\,Ho). 

The operator |l^|^/^i?o(^)|^|^''^ is trace class for all 2; G C off the positive real semiaxis 
(see [|5g, Problem 161]). This operator has limiting values \ V\^/'^Ro{E ± iO)\V\^/'^ for aU 
> in the Hilbert-Schmidt norm. Since the resolvent Ro{z) of Hq is an integral operator 
with kernel Ro{x, y; z) = ^^e*^!^"^!, it follows that lm\V\^/^RoiE + iO)\V\'^/^ has the 
integral kernel 

^\V{x)\'/\os{VE{x-y))\V{y)\'/' 

= ^^|V(x)|^/2 cosi^/Ex) cos{VEy)\V{y)\^/^ 
2v E 

+ -^\V{x)\'/'sm{VEx)smiVEy)\Viy)\'/', (2.11) 
2v E 
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and thus has a rank 2 and is obviously positive semidefinite. 



For positive energies the spectral shift function can be calculated as follows [|4^, [2^ ], 

(2.12) 



^(E-H ±\V\H)- ^ i„ det {l ±\V\^I-R,{E + 
m Ho ± \VIH,) - — log det(/±|V^|i/2i?„(i^;_,o)|V|V2) 



where the branch of the logarithm is fixed by the condition ^{E) ^ as — > cx). 
Obviously, the operator \V\^/'^Rq{E ± iO)|V|^/^ has no real eigenvalues (otherwise for 
some value of the coupling constant a G M the number —1 would be an eigenvalue of 
a| V"|^/2/2o(^±^0) I which implies that E e a^{HQ + a\V\) = 0, the positive singu- 



lar part of the spectrum of Hq + a\V\). Therefore, estimating the r.h.s. of ( [2.12D as in [ ]59| ] 
(Lemmas 2.2, 2.3, and proof of Theorem 3.1) for all > we obtain 

UE- Ho ± \V\,Ho)\ < c\\\lm\V\'/^Ro{E + iO)\V\'/^\\j, 



+ \\\V\^/^Ro{E + iO)\V\^/^\\^- 



J2 



where the constant C is independent of V and E. Here J7i and J'2 denote the trace and 
Hilbert-Schmidt norms respectively. 

Since lm.\V\^^'^Ro{E + iO)|V|^/^ > is positive semidefinite it follows that 

\\lm\V\^^^Ro{E + iO)\V\^^^\\j, = tTlm\V\^/^Ro{E + tO)\V\^^^ 

^ \V{x)\dx. 



l^fE 

Obviously, the r.h.s. is also a bound for WV^I'^Ro{E + iO)| 
Now we estimate ^ (0). According to ( |2.5| ) and ( |2.6| ) 



|e(0)|<no(i/) + -, 

where no{H) is a total number of eigenvalues of H. By the well-known Bargman-type 
estimate (see e.g. [^) 

no{H) < 1 + / \x\\V{x)\dx 



we obtain 



\m\<l+ [ \x\\v{x)\dx 



for E = 0. Since ^{E) is a nonincreasing function of < and ^{E) 
inf cr{H), the estimate ( PUD is valid for all E <0. 

Now let us fix some Eq > 0. For all E ^ [0, Eq] the estimate 



(2.13) 

for E < 



\m\<c- 



(1) 



(2.14) 
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with 



a 



(1) 



max < — h / |x| |V^(x)|(ix, 



C 



\V[x)\dx^ 



C 



\V{x)\dx 



is valid. The function ^{E) is continuous on [0, Eq] and thus by the Weierstrass theorem 
attains its maximum and its minimum, which by ( |2.10| ) and ( |2.13| ) are finite. Therefore, 
there exists Cy^ such that |^(-E')| < Cy^ for all E G [Oj-Eq]- This inequality combined 



with ( |2.14j ) gives the estimate ( |2.9| ) and completes the proof of the lemma. 



Remarks: 1 . The first term in the estimate ( [2.10[ ) with C = 1 /2 represents the high- 
energy asymptotics of the spectral shift function (see e.g. [ [TSi p3|]). Below (Lemma 6.1) 
we prove that for bounded potentials V with compact support the second term in ( [2.1U| ) can 
be omitted (with C independent of V). 

2. The estimate (|Z^) for > also follows from the results in [I 



Lemma 2.2 Let V be piecewise continuously dijferentiable, satisfy ( |2.i[ j and 

dV{x) 



{l + \x\) 

Then for each closed interval A C M 



dx 



dx < oo. 



where 



[ Ed^{E) = tT{V^/^E{A)\V\^/^ 

J A 



\ \ 1 dV(x) 

V{x) = V{x) + -X- ^ ^ 



(2.15) 



2 dx 

and E(-) is the spectral resolution for the operator H. 

Remark: A similar formula for the case of operators acting in LF'iW^) with d > 2 and 
A C (0, +oo) was proved by Robert and Tamura [^Tj] and by Jensen [^5p. 



Proof. Our proof closely follows the ideas of pT|]. Let /±(x, E) be the solutions of the 
integral equations 



f+{x,E)=e 



'wfEx 



sin V E{x — y) 



V{y)f+{y,E)dy, 



f.{x,E) = e-^^+ I ^^^^(--y) 



Viy)f4y,E)dy. 



These functions are solutions of the Schrodinger equation 

d^ 



dx'^ 



+ V{x)-E]fi{x,E)=0. 



(2.16) 
(2.17) 
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Let us also consider the functions 

,Pi{x,E) = T{E)f^{x,E), 

where T{E) is the transmission amplitude. They satisfy the Lippmann-Schwinger equa- 
tions which are integral equations of the Fredholm type ffl 



^P±ix,E)=^|J'^ 



(0) 



2y/E 



i\/E\x—y\ 



Viy)My,E)dy 



with 



are continuously differentiable with respect to cr. From ( |2.18| ) and since 



d 



-(a-V(x/a))^^^ = -2F( 



X) 



it follows that satisfy the integral equation 



d_ 



ij±{x,E]a) 



2VE 



dy. 



This easily gives 
d_ 



ij±{x,E;a) 



cr=l 



2 / R{E ±tO){x,y)V{y)ij^iy,E)dy, 



where R{z) is the resolvent of —d^/dx"^ + V . 



It is known [4-5] that the scattering matrix (2.7) can be calculated as 



S{E) = I 



2VE 



<i/{x,E)dx, 



where 



¥'>\x,E) = [i:f\x,E),ij^^\x,E) 
/ is the 2 by 2 unit matrix and * denotes Hermitian conjugation, such that 



(2.18) 



The functions ip± {x, E) are continuously differentiable with respect to x G M and E >Q 
[|18|,|15|]. Thus, 

E] a) = ip±{x/a, Ea'^) 



(2.19) 



(2.20) 
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Now we calculate 

d_ 



S{Ea')l=, = ^ [ (^a-'Vix/a)) ^^'\x, Ey^ix, E)dx 

2VE Jr \dcx J 

[ Vix)¥'\x,Ey^{x,E;a)U=idx. 



2VE 



By ( I2.19D and by the identity 

M^,E) = ij^^\x,E)- [ R{E ±iQ){x,y)V{y)^t\y.E)dy 

Jr 

we obtain 

-^S{Ea^) = 1= f V{x)'^{x,Ey '^{x,E)dx. 

OCT a=l sjE Jr 

Using the transformation property of the scattering matrix [ P3| ] 



^Ij+{x,E) \ _ ^^^^ ( ip-{x,E) 



we get 



± 

By the Birman-Krein theorem 



1 / dS{E)\ 



and therefore 

i\E) = E / nx)\Mx, E)\'dx 

for all E>0. 

Now we use the spectral representation for the spectral decomposition of the operator 
-d'^/dx^ + V (see e.g. \^), 

E{x,y,A)= ijj{x)ipj{y) 

■dE 
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where A+ = A n (0, +oo), A_ = A fi (— oo, 0) and tpj{x) are the eigenvectors of H with 
eigenvalues Ej < 0. 

Since V^/^E{A)\V\^/^ is trace class we have 



tr (vi/'E(A)|\/|i/2j = J2 I V{x)\i),{x)\Mx- I E^'{E)dE. 
It remains to prove that 

J2 [ VixM,{x)\'dx = - f Ed^E), 



j-.Ej&A 

or equivalently 



ti (V'/^E{A_)\V\'/^'^ = - f Ed^{E)= J2 (2-21) 



j-.EjGA^ 



The proof of ( |2.21D is quite elementary. It suffices to consider the case when A_ contains 
a single eigenvalue Eq of H with eigenfunction 00 such that E(A_) is the projector onto 
00- Let D be the generator of dilations, 



D = {2iy^ (x4- + -^x] . 

\ dx dx J 



In [27] it is shown that 



V = --[H,D]+H 

as a bounded operator from 11^^'^ (M) to 11^^'"^ (M), such that (0o, V4>o) is well defined. 
Then 

tr(i7i/^E(A_)|i/|i/2) = (0o,i?0o) =Eo- ^(00, [H,D]ct>o) = Eo. 



Lemma 2.3 For every z with Imz > the function logT(2;) is analytic and given by 

f i{E)dE 

\ogT[z) = - / — . (2.22) 

Jr rL — z 

Moreover, for all real X > Owe have 

log|T(A)|= / \og\E-\mE), (2.23) 
Jr 

where the integral is understood in the Riemann-Stieltjes sense. 
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Proof. The analyticity of logT(z) in the open upper complex half-plane C+ is well- 
known (see e.g. [jTSj]). Also logT(z) is bounded in C+ and 

hm logT{E + ie) = log \T{E) \ - 

e— >+0 

in all points of continuity of ^(E). Therefore, we can apply the Schwarz integral formula 
for the half-plane (see e.g. to reconstruct logT(z) from the limiting values of its 

imaginary part, 

logT(.) = -/i^ + .C, (2.24) 

Js. ^ — z 

with C being a real constant. For 2; ^ oo we have T{z) = 1 + / V{x)dx + 0(|2;|~^) 
Therefore, C in must be zero. ■ 



Below we will make use of the Aktosun factorizarion formula [|T]], which we formulate 
in the following form. 

Let V{x) be some real-valued locally integrable bounded function. Let {yn}n£Z be 
a sequence of real numbers such that y„ ±00 as n ^ ±00 and y„ < yn+i for all 
n G Z. Let Xnix) be the characteristic function of the interval [y„, Z/n+i]- We denote 
Vnix) = V{x)xnix) such that 



jyX) 



j=-n 



tends to V{x) as m,n — > 00. Let "'"'^ and Hj denote the Hamiltonians with domains 
of definition being the Sobolev space iy(^'^)(M), 

Let and Sj{E) be the corresponding S -matrices. 



We also consider the matrices 



T,{E) R,{E) 
L,{E) T,{E) 



1 R(-",rn) 



and 

A,(i?) 



1 

r,{£;) r,(£;) 

T,{E) T,{EY 
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From the unitarity of the scattering matrix (see e.g. [jTSJl) it follows that the matrices 
A^~"''''^\E) and Aj(E) are unimodular. We note that in Faddeev's terminology [ jTSl ] the 
elements of A{E) are given by the coefficients Cij{\fE). More precisely, 

C12(V^), -^ = -C22(v^), 



The Aktosun factorization formula states that 

m 

Ai-^'"^)(E) = Y[ Aj{E). (2.25) 

j=-n 

Here and below we understand the product Y[ in the ordered sense, i.e. 

n Ai(^) = A_„(E)---A^(E). 



The theorem below provides an alternative proof of ( [2.25D . Actually we show that the 



factorization property of the matrices A is directly related to the propagator property of the 
fundamental solution of the corresponding Schrodinger equation. 

Theorem 2.4 For arbitrary E > consider the matrix U{x, x'; E) G SL(2; C) which 
solves the initial value problem 

dU{x,x'-E) _ I y^^^M{x,E)U{x,x';E), U{x,x;E) = I, (2.26) 



dx l\fE 

with 

(1 p-2iy/Ex 
_e^^^^ -1 

The matrices A^^'^''^\E) are related to U (x, x'; E) such that 

A(-"'-)(E) = f/(y_„,y^;E). (2.27) 



The factorization formula ( |2.25| ) follows immediately from ( |2.27| ) and from the propa- 
gator property of U, U {x, x"; E)U{x", x'; E) = U {x, x'\ E). 

Remarks: 1. The matrix U (x, x'\ E) is related to the fundamental solution 0(x, x'\ E) 
of the Schrodinger equation 

d'^u 

+ V{x)u -Eu = 0, (2.28) 



dx"^ 
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which satisfies 



1 
V{x)-E 



(x, x'; E), x; E) = I. 



It is easy to see that 



where 



{x, x'; E) = P{x, E)U{x, x'; E)P{x', E)-\ 



P{x,E) 



-i\fEx 



i\fEx iy^e^*^^^' 



iV Ee 



such that 



1 

V{x) -E 



"^^^f'^^ Pfx, E)-^ + -^V(x)P{x, E)M(x, E)P{x, E)-\ 

ax r.- / 



2iVE 



2. The formula ( [2 .271 ) reduces the problem of the study of 4>{y-n, Vm', E) to the study of 
the scattering matrix for the corresponding single-site potential. Note that P(x, E)*P{x, E) 
is not a multiple of the identity operator. Therefore, the Hilbert-Schmidt norms of A*^""''") {E) 
and of (p{y-n, Vm] E) are not equal. This will become relevant in Section 5. 

3. The connection between the solutions of ( |2.26| ) and scattering characteristics was 
noted earlier in P7| ]. 



Proof. To prove the theorem it suffices to consider V{x) supported on the interval [0, a] 
and to show that U{x, E) := U {x, x'; E)\xi=q, i.e. the solution of the equation 



dU{x,E) 
dx 



2v^ 



V{x)M{x,E)U{x,E), U{0,E) 



satisfies 



R{E) 



U{a,E)=AiE)={ Ijlj 



T{E) T{E)* 



where T{E), R{E), and L{E) correspond to the potential V. 

Consider the solutions E) of the equation ( |2.28D , such that 



If jp\ - J e*^^' + i?(^)e-*^^, X < 0, 



^.ix,E) 



T{E)e 



isfEx 

1 

+ L{E)e 



\\fEx 



X < 0, 
x> a. 
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We introduce the 2 by 2 matrix 

W{x,E) = 
where A±{x, E) and E) are given by 

ip±{x,E) + 
ip±{x,E) - 



A^{x,E) = ^e'^^ 



A_{x,E) A+{x,E) 
B4x,E) B+{x,E) J ' 



i dip±{x,E) 



B^ix,E) = -e-^^^ 



\fE dx 

i dip±{x,E) 
\fE dx 



It is easy to see that W{x, E) satisfies the equation 

dW{x,E) i 



dx 



2VE 



V{x)M{x,E)W{x,E), 



(2.29) 



and 



Also we have 



W{a,E) 







L{E) T{E) 

Obviously, W{x, E)W{0, E)-^ also satisfies and equals / for x = 0. Thus, 

U{x,E) = W{x,E)W{0,Ey\ 

Moreover, 



U{a,E) 



since 



1 ^ \ f T{E) R{E) 
L{E) T{E) )\ Q 1 



T{E) T{E)* 

Finally note that U{x, x'; E), defined by (|Z2|), is given as U{x, x'; E) = U{x, E)U{x', E)'^. 



3 Cluster Property of the Spectral Shift Function 

In this Section we establish a cluster property of the spectral shift function for Schrodinger 
operators in (M) of the form 



H{d) 



dx'^ 



+ V,, V, = V, + V2{--d), 



(3.1) 
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where the Vi are in and have compact supports (see also ^ for results concern- 
ing cluster properties in the higher dimensional case). We study the behavior of the spectral 
shift function ^(E"; d) for the pair (H(d), Hq) when \d\ is sufficiently large. More precisely 
let V = V^Vi, V2) C M be such that the intersection of the minimal closed intervals con- 
taining suppVi and suppV2(- — d) is at most a point. This implies that Vi{x)V2{x — d) = 
a.e. for alld G V. We will henceforth assume that d E V. Denote 

CuiE; d) = i{E- H{d), Ho) - ^E; H^, H^) - ^E; H^, Ho) (3.2) 

with Hi = Ho + Vi (i = 1, 2). Also we set H2{d) = Ho + V2{- - d) such that H2 = H2{d = 
0). By the translation invariance of the spectral shift function, we have ^{E; H2{d), Ho) = 
^{E; H2, Ho) for aWd. ForbrevityinwhatfoUows we will write ^(E; d) = C,{E; H{d), Hq), 
UE)=mH,,Ho),t = l,2. 

Below we will need the following simple result: 

Lemma 3.1 Suppose that H^^\ /7q^' and H''^\ H^^^ are semi-bounded self-adjoint 
operators in the Hilbert spaces 7-^*^^^ and 7Y*^^^ respectively, such that H^'^^ — Ho \ i = 1,2 
are trace class. Then for a.e. E eM 

mH^'^®H^'\4'^®4'^) = 

^{E; © H^^\h^'^ © H^) + ^{E; H^^^ © H^^\ H^'^ © H^^^). (3.3) 

Moreover, 

^E- ff« © //f , h!,'^ © //f ) = m H^'\ h!,'^), 

aE; © if« © iff) = m H^'\ H^o^) (3.4) 



a.e. on 



Proof Let Ti = H^^'^ © H^^l For every / G C^{M) we have 
tr^ ©//(^))-/(if«© iff)) 

= tin [fiH^'^) © f{H^'^) - fin!,'^) © fi4'^) 
= tr^(i) {f{H^'^) - f{Hl>'^)) + tr^(.) {f{H^'^) " f{H^o 
= tin {f{H^'^) © /(/if) - /(iif ) © fiH, 

+trn (/(iif ) © fiH^'^) - /(iif ) © fiH, 
= tr^(/(i7W©iff)-/(iff ©iff; 
+tin (/(i/f ©//(^))-/(i/f ©ifi 



(2)^ 



r(2)- 
'0 , 



From this it follows that ( |3.3| ) and ( |3^ are valid up to an additive constant. From our 



normalization and from the semiboundedness of the operators involved it follows that this 
constant equals zero. ■ 
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Theorem 3.2 For all d e V{Vi, V2) 



Proof. Let us fix some y E M. lying between the supports of Vi and V2(- — d). More 
precisely we require that y is such that suppV^i C {—00, y] and suppV2(- — d) C [y, 00) 
if > and suppV2(- — d) C {—oo,y] and suppVi C [y, 00) if c? < 0. Without loss of 
generality we may suppose that d > 0. For every potential V satisfying ( p?T| ) along with 
the Hamiltonian H we consider the self-adjoint Schrodinger operators on L^(R) 

with Dirichlet and Neumann boundary conditions atx = y.ln accordance with the decom- 
position 

L^R) = L\-oo, y) © L^{y, 00) 

one can write if (^•^) = jj^'^^ H^'^K From Krein's formula (see e.g. [Il4|]) it follows 
that H^^^ and H^^^ are rank one perturbations of H. Also we have 

if(^) <H< /J(^) (3.5) 

in the sense of quadratic forms. Since Vi and V2(- — d) have disjoint compact supports, we 
have 

From Lemma 3.1 it follows that 

= m '"^ i^i"'"^) + Hr\d), (3.6) 
Now by the chain rule for the spectral shift function and from ( |3.6| ) it follows that 



= m Hid), iJ(^'^)(rf)) + i^r''^^ i/o) 

iff'^^ if^'"^) + m i^r'"^ (3.7) 

On the other hand and again by the chain rule one has 

+m i^r'^\ i^i^'^^) + m i^r-'^^ ^0) (3.8) 
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and analogously with Hi replaced by H2{d). From ( P77| ) and ( |J^ ) it follows that 

^E; H{d), Ho) = ^E; H^, Ho) + ^E; H,{d), Ho) 
H{E; H{d), - m Ho) 

-m H,, i/f'^)) - m H^id), iff 

such that 

ei2(i?; d) = -^E; H{d)) - m ^r'^^ ^o) 

H{E; Hf'''\ H,) + ^E; iff '^^(rf), H,{d)). (3.9) 

We note that in the terminology of Gesztesy and Simon [gl]] ^{E; H^^\d), H{d)), 

^{E; hI^\ Hi), i = 1,2, and ^{E; Ho^\ Hq) are the xi-functions for the operators H{d), 

Hi, i = 1,2, and Ho respectively. 

From (P3|) and from the fact that the absolute value of the spectral shift function for 



rank one perturbations is not greater than one, we have that for all real E 

< ^{E;H^^\H) < 1, 
-1 < ^{E;H^^\H) < 0, 

where H stands for one of the operators H{d), Hi, H2{d), or Ho- Also (see pl|]) 

.(p.H^D) rr. i 1/2, >0 

aE,Ho - I , 

aE,Ho ,^o) - I . 

Therefore, from ( p.9D we obtain 

-3/2 <^i2{E;d) <3/2, E>0, 
-1 <^i2{E;d) <2, E<0 

using Dirichlet boundary conditions and 

-3/2 < ^i2(^;rf) < 3/2, E>0, 
-2<ii2{E-d) <1, E<0 

using Neumann boundary conditions. This completes the proof of the theorem. ■ 

For our further purposes to be elaborated below Theorem 3.2 provides a completely 
sufficient information. However, we note that when E > the bound for ^12 (-E; d) can be 
improved: 

Theorem 3.3 For all d E V{Vi,V2) and all E > 

\^i2iE-d)\<l/2. (3.10) 
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The proof of Theorem 3.3 will given below. 

Theorems 3.2 and 3.3 imply that \^u{E; d)\ < 1 for all d G V{Vi, V2) and all E e R. 
Let us define the functions 



^^\e) = ^,iE)±l,j = l,2. 



(3.11) 



Corollary 3.4 For all E the functions ( \3.11\ ) satisfy the inequalities 

ly{E-d)<ip{E)+ip{E), 
e-\E;d)>^[-\E)+^i-'\E), 



(3.12) 
(3.13) 



i.e. ^^~^\E) is a subadditive and KE) is a superadditive function with respect to the 
potentials Vi and V2{- — d). 



Proof. Due to (|3.2[ ) and Theorems 3.2, 3.3 we have 

HE- d) = UE) + UE) + ii2{E- d) < UE) + UE) + 1. 
Adding 1 to both sides of this inequality we arrive at ( |3.12D . Similarly, we have 

i{E- d) = UE) + UE) + ii2{E- d) > UE) + UE) - 1. 
Subtracting 1 from both sides yields ( |3.13D . ■ 



Proof of Theorem 3.3. Recall (see [[T^, |T5]1) that the scattering matrix at energy > 
for the pair of Hamiltonians {Hj = Hq + Vj, Hq) is given by 

c.(j^) _ f Tj{E) Rj{E) \ ■ y 



We use the fact that for all d e V{Vi, V2) and all E > 

^i2{E;d)-- 



1 , 1 - Ri(E)L2(E)e^'^'^ 
log 



27ii ''I - R^{EYL2{EYe-'^^^'^' 



(3.14) 



The formula ( p.l4[ ) has appeared earlier in [ |5^ and follows from the Aktosun factorization 
formula [2.25| 



T{E-d) 



T^{E)T2{E) 



I - Ri{E)L2{E)e^i^<^ 



(3.15) 
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and the identity (see e.g. p5|]) 



1^ = det S,{E) = e-^-^^(^), J = 1, 2. (3.16) 



According to ( |2.7D the reflection amplitudes Lj{E) and Rj{E) can be represented in 
the form 

(E) = A, (E) e'^f^ , Rj (E) = A, (E) e'^f^ 

with < Aj{E) < 1. Moreover Aj{E) = 1 only when Tj{E) = 0. This in turn can 
happen only if E = (see [[TEI, [T3|]). Therefore 



1 - R^{E)L2{E)e^'^'' 



log 



1 - Ai(E)A2(^)e-*(4'''+4^^+2v^'i) 



-2i arctan ■ 



1 - Ai(E)A2(^) cos(5f ^ + + 2VEd) 

Due to the fact that < Aj (E) < 1 for > the absolute value of this expression is 
strictly bounded by vr. ■ 

Since the spectral shift function is right continuous the estimate 

\^i2{E-d)\ < 1/2 (3.17) 



holds also for E = 0. We note that the inequality ( p.lVD for = also follows from 



the Levinson theorem |^], Theorem 3.1 of p3| ] and the results of [0]. We include this 
discussion (see also [|54|]) to show the different aspects of the problem. For brevity we set 
n(y) = noiy), the number of bound states below E = 0. Theorem 3.1 of [ |33| ] states in 
particular that n{Vi + V2(- - d)) = n{Vi) + n(V2(- - d)) = n{Vi) + n(y2) for aU d such 
that d E ViVi, V2) if = is an exceptional point for at least one of Hi or H2. If = 
is a regular point for both Hi and H2 then 

niVi) + niV2) - 1 < n{Vi + V2(- - d)) < niVi) + 72(^2)- 
Therefore, by ( |2.5p and ( P^ we have 

-ei(0) - 6(0) < n{Vi + V2{- - d)) < -6(0) - 6(0) + 1 
if £■ = is a regular point for both Hi and H2, and 

n{Vi + V2{- - d)) = -6(0) - 6(0) + 1/2 



if _E = is an exceptional point for at least one of the operators Hi and H2. Then according 
to Theorem 2.3 of [@] we have the alternatives 
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(i) Let = be an exceptional point for exactly one of the operators Hi and H2. Then 
E = Ois a regular point for H(d) and therefore 

e(0;d)=6(0)+6(0). 

(ii) Let £■ = be an exceptional point for both operators Hi and H2. Then = is an 
exceptional point for H(d) and therefore 

e(0;d) = ei(0)+6(0)-l/2. 

Let now _E = be a regular point for both Hi and H2. The discussion at the end 
of Section 2 in [0] states that = is a regular point for H(d) for all except possibly 
one value d = do e R . Excluding this value do (however, we found no proof that do ^ 
V{Vi,V2))we get that 

ei(o) + 6(0) - ^ < m d) < 6(0) + 6(0) + \. 
for all d e V{Vi, V2) with d ^ do- 



Now we reconsider the case E < 0. We recall that in this case the spectral shift function 
equals minus the number of eigenvalues less than E. The discussion below shows that the 
bound |62(-E^; c?) | < 1 for < in general cannot be improved. 

Let i?^*^ < 0, ki = 1, . . . , no{Vi), i = 1, 2 be the eigenvalues of Hi. Let Ek{d) < 0, 
k = 1, . . . , no(Vi+V2(-— (i)) be theeigenvalues of if((i). We recall some known properties 
of Ek{d) g. The functions E^id) are real analytic ind e V = V2). First let 

< be an eigenvalue of both Hi and H2. Then H{d) has two eigenvalues E±{d), 
E^{d) < E < E+{d), which both converge to as 00. Moreover, E'_ > and 
E'^ < (with prime denoting the derivative with respect to d) for all d G V. Secondly, 
let E be an eigenvalue of (say) Hi but not of H2. Then the only eigenvalue E{d) of H{d) 
approaches E as d ^ 00. Moreover, either E'(d) > 0, or E'{d) < 0, or E'{d) = for all 
d & V. There are no eigenvalues of H(d) other than the ones described above. We note 
also that for all d G V either -Efc((i) 7^ E^^'' for any k, k^ and i, or Ek{d) = E^^^ for some 

k, ki and i. Indeed let us suppose that there is do G V such that Ek{do) = E^^^ = Eo for 
some k, ki and i. Then inspecting the proof of Theorem 1 .2 in P3| ] we see that this implies 
Ek{d) = Eo for all d E V. As is well known in the one-dimensional case the eigenvalues 
of H{d) are simple and therefore the Ek{dys cannot cross each other. 

Now fix some Eo < 0. There are two cases to be considered: (i) Eo is an eigenvalue 
of neither Hi nor H2, (ii) -Eo is an eigenvalue of at least one of the Hamiltonians Hi, H2. 
By the discussion above in case (i) there is d{Eo) > such that for all d > d{Eo) one has 
f^Eoi^i ~^ ^2{- — d)) = UEoiVi) +^£'0(^2). When d G V decreases only one of the curves 
Ek{d) can pass through Eo (since the Ek{dys cannot cross each other) thus decreasing or 
increasing nEo{Vi + V2(- — d)) by one. In case (ii) nE^iVi + V2{- — d)) does not depend 
ond eV and equals UEoiVi) + nEo(V'2) or differs from nEo{Vi) + nEo{V2) by one. 
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4 Existence of the Spectral Shift Density 

As stated in the Introduction we will consider random Schrodinger operators H{uj) in 
L^(M) of the form ( |1 . 1| ) with {aj{uj)}j^z being a sequence of i.i.d. variables on a prob- 
ability space JF, P) having a common density which is continuous and has support 
in the finite interval C M. Also the sequence {aj{uj)}ji^z is supposed to form 

a stationary, metrically transitive random field, i.e. there are measure preserving, ergodic 
transformations {Tj}j^z on Q such that aj{TkUj) = aj-k{u!). 

We suppose that the single-site potential / is in C(M) with supp/ C [—1/2, 1/2] and 
/>0. 

First we introduce the Hamiltonians 

= + E %(^)/(- -^0, (4.1) 

j=-n 

such that i/(")(c<;) = i7(-"'")(cj). Let u;) be the spectral shift function for the 

pair (a;), i/o)- 

The operators H^~^''^\TkUj) and ("""^'™"^)(ct;) are unitarily equivalent. In fact, con- 
sider the unitary shift operator f/j, j G Z on L^(]R) given as Ujf{x) = f{x — j). Then one 
has 

j=-n 



j=-n 
j2 m—k 



j=—n—k 

= [/fci/(-"-^'"-^)(cu)f/fc*. (4.2) 

Since the spectral shift functions for pairs of unitarily equivalent operators are equal, we 
have 

This remains true for the functions 

to) = ^i-n-k,n.-k)^^. ^) ^ 

Now let A; be an arbitrary integer such that —n < k < m. Then due to Corollary 3.4 
we have that 

ei""'"'^ {E;uj)< a;) + '"^ ^) (4.3) 

and 
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Now we show that 

r+ = inf J E\^i-'"'\E;uj)]>-oo, 

m,n m + n + 1 I J 

r- = sup J ^ E{t''^"'\E;u;)]<oo, 
m,n rri + n + l I J 

where E denotes the expectation with respect to the probability measure P. 
First we note that 



< 



From the inequality (4.3) it follows that 



uj)<J2 Ho + a,{uj)f, Ho) + {n + m + 1] 



]=-n 



and then by the monotonicity theorem of the spectral shift function with respect to pertur- 
bations we have 

et"''"^^; oj)<{n + m + l) UE- Ho + «+/, Ho) + 1] . 
Hence by Lemma 2. 1 

T^<i{E-Ho + a+f,Ho) + l< oo. 
Similarly we can prove that r+ > — oo. Indeed, 

r+ = inf^i--E{^l-"'-)(E;^)| 

m,n m + n + 1 L ) 

= inf ^ E|d'"'"^(E;^) + 2| 

m,n m + n + I I J 

> inf E|ei~"''"^(^;cu)| 

m,n m + n + 1 I J 

m 

> mf^——-Y,i^{^iE;Ho + a,{u),Ho)}-l) 
m,n m + n + 1 ^ — ^ 

j=-n 

> aE;Ho + a.f,Ho)-l> 
Thus we have proved 

Theorem 4.1 For every E e R the family ^^^'''"'\E; Lu) is a subadditive and "''"^ {E; lu) 
is a superadditive random process. 
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Applying now the Akcoglu-Krengel superadditive ergodic theorem we obtain that for 
every E eM. there is a set ^7^; C of full measure such that 

lim ^—-^^=:m (4.4) 

exists and is non-random. We call this limit the spectral shift density. 

Now the problem is to show that the set VIe of full measure can be chosen to be inde- 
pendent of E as long as is a point of continuity of the limit. We note that a priori the set 
n Em^E is not necessarily of full measure. 

We recall how this problem is solved for the density of states N{E) (see e.g. [ [T3| , p. 
312]). Once one has established the existence of the limit 

lim tl^ \f^ = N{E) (4.5) 

m,n^oo n + m + 1 

for every fixed E and P-almost all uj E Vt, one can choose Vt as the intersection of all sets 
VLe when E runs through the rationals and redefine the limiting function N{E) to make 
it right continuous. Since N{E) is a monotone nondecreasing function of E, this could 
change the values of the limiting function on at most a countable set of discontinuities. 
Hence ( |4.5P is valid at every continuity point of N{E) for all uj eVL, which is obviously of 



full measure. 

In our case the limiting function ^ {E) is not monotone. However intuition says that 
i{E) must be equal to Nq(E) — N{E) (this is indeed the case as will be proven below). 
Therefore, ^ (E) is expected to be at least of bounded variation. The simplest way to prove 
this is to show that {n + m + l)~^^i~'^''"'' (E) are Lipshitz functions with Lipshitz constants 
bounded uniformly in n and m, which is nothing but a Wegner-type estimate for the 
spectral shift density. This guarantees that ^(E) is Lipshitz continuous. 

We will need the spectral averaging theorem [I 



Lemma 4.2 Let H = Hq + W with W e L^{R). Let V e L^(R) be nonnegative. Es(-) 
the spectral decomposition of unity for Hg = H + sV. Then for any Borel set /S. <ZM. 

r ti{V'/^Es{A)V'/')ds = [ i{E-Hs,,Hs,)dE. 

J so J A 



Remark: ti{V^^'^Es{A)V^^'^) is well defined since for every / G L^(R) the operator 
/i/2^,(A)|/|V2 is trace class (see [0). 

The present formulation of the spectral averaging theorem is a direct consequence of 
a slightly extended version of Theorem 4 in [^. This extension is straightforward and 
therefore we do not discuss details here. Now we prove 
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Theorem 4.3 Let f be piecewise continuously dijferentiable with supp/ C [—1/2, 1/2] 
and assume there is a constant Cf > such that 



f{x) + 



X df 
2dx 



< Cffix) 



for a.e. x G supp/. Then E^{E) is Lipshitz continuous for all E eM., i.e. for each closed 
interval A C M there is a constant Ca such that 



\E2^{E2)-E^^{E^)\<Ca\E2-E, 



for all El, E2 e A. 



Thus it suffices to determine ^ (E) for E running over a Lebesgue-dense set of M, say 



the rationals. We choose the set Q as the intersection n^;^^^^^;. The limit then exists 
for all E G M and all u eU. 

Now we turn to the proof of Theorem 4.3. We note that the method of the proof is not 
restricted to the one-dimensional case and can be extended to higher dimensions (see [06|]). 



Proof of Theorem 4.3. We apply Lemma 2.2 to if = ii'(-"''")(a;). Let us set Ai 
[£"1, £2] ^ A. Then we have 

E2i^-^^^\E2;uj) - Eii^-^''^\Ei-uj) 

i^-''^"'\E-uj)dE + [ Edi^-''^'^\E-uj) 



Ai 

m 



n II tj 

^1 j=-n 

I ^(--'-\E;uj)dE - Y «o(T_,a;)tr 

"^^1 j=-n 



where 



and E 



(— n.m) , 



is the spectral resolution for H^~'^''^\uj). Therefore 



E { ^2^^-"'"^) (^2 ; ^) - Ei^^-"'™) {El, 00)} 

= e|^ f]E{aoMtr(7i/2Ei--"-^'-^^^^ 

First let us estimate the second term on the r.h.s. of (|4.6D, 

E {aoMtr fr/2E(7"-^-'-^-)(Ai)|7|V2 



(4.6) 



< a+E 



tr (/i/2E(7"-^'-^-)(Ai)|/r/') 
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Since for any A E J'l the inequality |trA| < tr|y4| holds, we have 

tr (7V2E(;"-^'™-^-)(Ai)|7|V2)| < tr (|7|i/2E(7"-^'-^-)(Ai)|7r/' 

<CftT (/V2E(;"-^'-^)(Ai)f/2). 



Let us denote 



j=~n 



and let Ei,,^ (■) be the corresponding resolution of the identity. Then we have 



E{tr (/V2E(- 



n~j,m-j) 



(Ai)0} 



E 



< 



Now we apply Lemma 4.2, according to which we obtain 

Pdatr(/V2E(;„"-^'™-^)(Ai)/^/') 

""-^■''"-^■)(cj)), 



where ^{E; Ha+ 



Ai 

(-n-i,m-j) 



/ dEaE;Hi~;^-^'^-^\u),Hi- 

JAi 



{-n~j,m-j) 



cj)) Stands for the spectral shift function of the 



pair {Haj^ ^'"^ ^\^), Ha^ ^\uj)). By the chain rule and Corollary 3.4 we have 



aE-H^---^^^-^\uj),H^^ 



-n~j,m-j). 



n-j,m-j). 



0) 



< m HtV'"'-'\uj), Ho) + m Ho + a^f, Ho) + 1 

-m Htl-''"'-'\^), Ho) - m Ho + «-/, Ho) + 1 
= ^{E-Ho + a+f,Ho + a.f) + 2. 

Therefore, the second term on the r.h.s. of ( ^^61 ) can be bounded by 

(n + m + 1)1^2 - ^il II0IIOO max ^{E; Ho + a+f, Ho + a.f) + 2 



(4.7) 



By Lemma 2.1 and since ^{E; Ho + a+f, Ho + a-f) = ^{E; Ho + a+f, Ho) - ^{E; Ho + 
a^f, Ho) the maximum of ^{E; Ho + a+f, Ho + a~f) is bounded. 

Now we estimate the first term on the r.h.s. of (P^). By the Fubini theorem 
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Using monotonicity and Corollary 3.4 we obtain 

E{^(-"''")(i?i;^)} < {n + m + l)[aE,;Ho + a+f,Ho) + l]. (4.8) 

The expressions in the square brackets in ( |4.7D and are finite by Lemma 2.1. 
Thus we have proved that 



(n 



+ m + 1)-' \E2E {e(-'^''"H^2; oo)} - E,E {^(-"'"^^^1; ^)}| 

<Ca|^2-^i|. (4.9) 



Now we note that by the Lebesgue dominated convergence theorem 

lim (n + m + l)-iE{ei""''"n^)} =^(^) 

m,?i— >oo 

for every fixed G M. Thus, taking the limit n, m — > 00 in (P3|) we arrive at the claim of 



theorem. ■ 

Theorem 4.4 The spectral shift density is the difference of the integrated densities of 
states for the free and the interacting Hamiltonians, C,{E) = Nq{E) — N{E) for all E G M. 

Proof Let D*^~"''"^(ti;) and Dq""'*"^ be the self-adjoint Schrodinger operators corre- 
sponding to the differential expression ( [4.1[ ) and Hq = /dx^ respectively on L'^(—n — 
1/2, m + 1/2) with Dirichlet boundary conditions at x = —n — 1/2 and x = m + 1/2 
(the Friedrichs extension of ( |4.1D on C^{—n — 1/2, m + 1/2)). They have purely discrete 
spectrum and therefore the spectral shift function ^{E; D^~"''"^\uj) , Z^q^"'™^) is simply the 
difference of the corresponding counting functions 

NiE; dJ""'""^) - N{E; D^"'^"'\uj)) . 

It is well known (see e.g. [[131]) that for all G M 

m,n-~>oo n + m + 1 
m,n->oo n + m + 1 

almost surely. Hence 

N,iE)^NiE)= lim «^;^—'(^).^^—') (4.10, 

m,n^oo n + m + 1 

for almost all G Vl. 

Now we prove that the difference 

C{E- H(--'"^\u;), Ho) - m D(--'"^\u;), D^-''")) 



4 Existence of the Spectral Shift Density 



28 



is bounded in absolute value by 2 uniformly in m G M and lu E il. Therefore, the 
existence of the limit ( |4.10| ) immediately implies the existence of ^(E) and the equality 
$,{E) = No{E) ~ N(E). Thus, in fact we do not need Theorem 4.3 in this context. 
By the chain rule for the spectral shift function we have 

H{E-,Hl,--'-\u),H;^-'-'>)+aE-,H^-^'"'\Ho). (4.11) 



Here H^^ "''"^(u;) denotes the operator ( |4.1D on with Dirichlet boundary conditions 



at X = —n — 1/2 and x = m + 1/2 such that 

/7(^-"''")(cu) = D^-"^'-"^ © © (4.12) 

with respect to the direct decomposition of the Hilbert space L^(]R) = oo, —n — 
1/2) ®L^{-n -l/2,m + 1/2) ©L^(m + 1/2, oo). Similarly there is a decomposition of 



H, 



(~n,m) 
0,D 



with the same boundary conditions. From Krein's formula it follows that "''"^\uj) and 
Hq~^'"^^ are rank two perturbations of i/*^^" ™) (cj) and Hq, respectively. Thus, we imme- 
diately have 

< ^{E;H^j^'''"'\uj),H^-"'^"'\uj)) < 2, 

0<aE;H^^j^'"'\Ho)<2 (4.13) 

for all n, m G N and cu E Vl. Actually, ( |4.13| ) can be improved [ pS] , Remark 5.2] such that 

l<mKD"'\Ho)<l. 



From ( |4.12| ) by Lemma 3.1 it follows that 

aE;H^-''''^\uj),H^-^-"''>) = e(^;Z^(-"''")M,Z}i""''")). (4.14) 
This completes the proof of the theorem. ■ 

Remark: We comment on the formula ( ^~U\) . We note that ^{E; iJ(-"''")(u;), Hq) is 
continuous in E > although C.{E; H^'^'"^\ Hq~^'"^^) is a step-like function being the 
difference of two counting functions. This means that the difference 

^E; Hi^^''^\ Ho)) - m Hi,--'-\u), H^~-'-\u)) 

compensates the jumps of C{E] H^^'"^\uj) , H^'j^'™'^) making the r.h.s. of ( |4.1 1| ) contin- 
uous. This was noted by Jensen and Kato [|2^]. A similar phenomenon was found for 
obstacle scattering in by Eckmann and Pillet [|T7|]. 
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Fig. 1: Spectral shift density for the deterministic Kronig-Penney model (see text). 

As an illustration to Theorem 4.4 (see Fig. 1) we have calculated (n+m+l)~^^'^~"''") (E) 
(dotted line) for n = m = 7 with / taken to be the point interaction and aj{uj) = 1 (the 
Kronig-Penney model) and compared this result with No{E) — N{E) (solid line). The 
density of states N{E) for this case can be given in closed analytic form (see e.g. [|]]). 

5 Density of the Transmission Coefficient 
and the Lyapunov Exponent 

Now we turn to a discussion of the density of the transmission coefficient. Let j^^""'™) 
be the transmission amplitude at energy E > for the Hamiltonian H'^~"'''^\uj) ([0]). We 
now prove 

Theorem 5.1 For every fixed E > and almost all uj eVL the limit 

n,m^oo n + m + 1 
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exists and is non-random. 

Remark: For periodic deterministic potentials the behavior of T^""'"^ {E) as n ^ oo 
was studied numerically in [53p ■ 



A very similar statement was proven earlier by Marchenko and Pastur p20 . Our proof 



is a slight modification of that given in [|42| 1 . We start with 



Lemma 5.2 For all E > the transmission amplitude for the Hamiltonian ( |3J| j satis- 
fies the inequality 

\T{E)\ > ^-\T,{E)\ \T,{E)\. 



Proof. By the Aktosun factorization formula 

T{E) = T^{E)T2{E) (1 - R,{E)L2{E)r' . 

By the unitarity of the scattering matrix \Rj(E)\ < 1 and \Lj(E)\ < 1, j = 1,2 for all 
E>0. Hence 



\1 - Ri{E)L2{E)\ < 2 

and the claim follows. ■ 

Proof of Theorem 5.1 We set 

t^^^\E) = Uog\T^--^^\E)\<Q. 

From the fact that 



which is is an immediate consequence of ( |4.2[ ) we obtain 
By Lemma 5.2 we have 

tl-'''''\E) + t^^^'^\E) < t(7"'™)(^) 

for all k with —n < k < m. Thus t\^^'^\E) is a superadditive random process. Theorem 
5.1 now follows by Akcoglu-Krengel superadditive ergodic theorem. ■ 

Theorem 5.3 For every E > 7t(-E) equals 7(-E), the upper Lyapunov exponent for 
the fundamental matrix of the Schrddinger operator 
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This observation (although without a complete proof) is known (see ^]). We start 
with recalling the definition of the Lyapunov exponent (see e.g. p8|]). Let (f)i^{x] E) be the 
fundamental matrix of the Schrodinger equation 



such that 



+ Y,(^Mf{--j)i^ = Ei^ (5-1) 



for any solution of (|57T|). The Lyapunov exponents 7^(-E) are defined by 

7^(i?)= ^log||</..(x;i?)||. (5.2) 

Actually it can be shown that 'y^{E) = 7j(-E) = 'j{E) for fixed E and P-almost all u. 
Moreover, lim can be replaced by hm (see e.g. p8|]). 

Now fix E > 0. We have to show that for P-almost all tu 

7(E) = lim ^log||0,(x;E)|| =7t(^). 

First let us redefine the fundamental matrix such that 

u'{x + l/2) J '^-^^'^^ 1^ u'{l/2) 

This obviously does not change the Lyapunov exponent, which we can calculate as a limit 
over integer x, 

7(E) = Urn -logU^{n;E)\\. 

This implies that we can express (puiin; E) in terms of the transmission and reflection ampli- 
tudes for the Hamiltonian //(i'") (u;) . More precisely, let us consider the particular solutions 
il)^^\x] E) of the Schrodinger equation with Hamiltonian H^^'"'\tu) corresponding to the 
energy E > and for x < 1/2 having the form 



^W(a;;E) 



e 



Then it is easy to see (see [ [TSl ] and Section 2), that for x > n + 1/2 these solutions have 
the form 

lL ' -* (E) ^ 1 ^-i^x 



i^^J-\x-E) = ^pL-\x;E), 
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where t!:!''^\e) and L'^''^\e) are transmission and reflection amplitudes for t 
nian H^^'"\uj), respectively. Therefore, the matrix elements of (p^jin; E) are g: 



1 
2 



1 

+ 2 



1 



6 ~r 7"; ^ I 



12 



1 



I/j' ''(-5')* _ ^^i/_^C^giv^(n+l) 



21 



2 



:^^il^(^Tg-iVS(n+l) _ (-^) ^iVg(n+l) 



0a; (n; -E) 



22 



1 
"2 

1 

2 



^i\/En ^—i\/~En 

W\ey~tF\e) 

_g-jv^(n+l) _^ iV^(n+l) 



Now using the relation (see 



of 4)^{n;E). 



we calculate the Hilbert-Schmidt norm 
we get that 

i;E)\\% = \Ti''-\E)\-'{2il-\Ti''-\E)\')cos'e, + 2cos'9, 



in: 



After some simple transformations 



+E 



|Ti''")(E)|2 sin 01 - sinks' 
1- |Ti''")(^)|2sin0i + sin02^ }. (5.3) 



Here for brevity we have introduced the notations 

e, = v^(n + l)-0(i'")(E) + 7r/2, 

Obviously, the expression in the braces on the r.h.s. of ( |5.3D is bounded from above for 
every E > uniformly in G N by 4(1 + E + 1/ E). In Appendix A we show that this 
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expression is also bounded from below by the positive constant 

AE 

uniformly inn E N. Thus, 

7(E) = lim -log||0^(n;E)|| 

n^oo n 



Now consider 

m,n) 



- lim i|og|r^">(£;)|. (5.4) 



n,m^oo n + m + 1 

iog|ri^"'='^''=)(i?)| 

— lim lim 



, , log 14^(^)1 
= — lim iim , 

where N = n + m + 1 and k = 1 + m. By for P-almost all u 



independently of k. 



Now we recall that argT^ "'™)(E) = -7r^(~"''^)(E; w). Then by (g]|) and Theorems 
5.1, 5.3 we have 

Corollary 5.4 For every E > and F-almost all to 

lim i^^^^^— if^ = -7(i?)-.<(i?). (5.5) 

m.,n->oo m + n + 1 



We note that Corollary 5.4 can be reformulated in such a way that it permits a general- 
ization to the higher-dimensional case. We recall that (see [^3]]) 

Ti-"'™)(E)-i = det (/ + Vj-"'™)^^^i?o(^ + ^0)|rj-"'™V^^) , (5-6) 

where 

m 
j=-n 
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Now we can rewrite (p3|) as follows 



lim 



logdet (^/ + \/j-"''^)^^^/?o(^ + iO)|yj-"''"V^^ 

= 7(^) + ^<(^). 



In this form ( |5.5D can now be generalized to the higher-dimensional case, thus defining 
'y{E), which is something like the multi-dimensional Lyapunov exponent (for details see 
[B^). Also note that the above determinant equals the determinant of the lost matrix [B5p. 



We cannot use ( p3| ) to calculate 7(-E') for < 0. However, this can be done by means 
of analytic continuation (see Section 6). 



We turn to the claim '-/{E) > for almost all E > 0. We start with some preparations. 
Let us denote 



By the identity 



1 












Li-"'"'\E) 




1 


Tf'""\E) 




"•'"'(E)* 



00) 



|2 



1^2 |Ti-"'™)(E)|2 |Ti-"''")(E)|^ 



and by Theorem 5.3 one has that for every E > Q 



l{E)= lim / log||A(-"'-)(E;cu)|| (5.7) 



almost surely. Also ( |5.7| ) follows directly from the definition of the Lyapunov exponent 
(O) and Theorem 2.4. 



Let i^Q, := Hq + af for some a e M. The corresponding elements of the scattering 
matrix at energy > we denote by Ta{E), Ra{E), and La{E). Let 



AaiE) 

and 



1 RcjE) 

Ta(E) Tc.{E) 
Lc{E) 1 

Tc,(E) Tc,{E)* 



-e 



Obviously, 



A,(E;a;) = [/^A,^H(E)t/^^ 
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where 



i^/E 



Up 







e 

From the Aktosun factorization formula ( [2.25P it follows that 

m m 



j=-n 



■m-1/2 



(5.8) 



Since t/g is unitary one obtains 

||A(-"''")(E;u;)|| = 

with 



n ^-.h(^) 



rc(E) 



' T^(E) 



Tc{E) T^(E)* 



(5.9) 



Since the aj{uj) form a sequence of random i.i.d. variables, for every E > the sequence 
Aajiuj)iE) is a sequence of random i.i.d. SL(2; C)-valued variables with corresponding 
distribution (pE- Recall that the distribution density is supposed to be continuous with 
compact support. 

The matrices Aa{E) have the form 

a b 
b* a* 

The closed subgroup of all matrices from SL(2; C) having this form (which we denote by 
SL]k(2; C)) is isomorphic to SL(2; M). Indeed, 



a b 

b* a* 



Q 



Rea + Re6 —Ima + Imb 
Ima + Im6 Rea — Re6 



Q 



-1 



with 



^ 2 



1-i 1+i 
I- I -I - I 



Let ia{E) be the spectral shift function for the pair of Hamiltonians {Ha, Hq) such 
that TaiE) = |T„(E)|e-*^«"(^) for E > 0. It is well known (see e.g. 0) that r^(^), 
Ra{E), and La{E) at fixed energy E > are real analytic functions of a G M. Since 
/ has compact support, Ta{E), Ra{E), and La{E) at fixed a G M are real analytic with 
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respect to E > [jT5|]. Moreover, they are jointly real analytic in a G M and E > 0. Since 
tS{E) for all E > 0, |T^(^)| and are also jointly real analytic in a G M and 

E>0. 

We recall that potentials with compact support cannot be reflectionless, i.e. Ra{E) = 
for all E > implies a = [|l5|]. 
By real analyticity the set 

S^ = {E>0: R^{E) = 0} 

for every fixed a 7^ is discrete or empty. Let 

S = {E>0: Ro,{E) = for all a G suppy?} . 

By the assumption that supp(^ has a positive Lebesgue measure and by real analyticity of 
Ra{E), the condition Ra{E) = for all a G suppy^ or even for a in a subset of suppv? 
of positive measure implies that Ra{E) = for all a G M. Obviously, S = CiamSa and 
therefore is also discrete or even empty. 



Remark: One can easily show that a necessary (but not sufficient) condition forEeS 

is 

We know, however, no example of a potential with 5 7^ and expect that actually 5 = 0, 
since, intuitively, it is clear that for the potential / > at hand Ta{E) — »• as a — 00. 



Now for every fixed E > 0, E ^ S we define the functions 



rJe)^^^) 

Let 

5'(±) = |e > 0, E ^ S : F^\a) does not depend on a G suppv^j . 
If for some E > one of the functions = C*^^^ = const for all a G supp(y9, then 



cos2(/E-<„(E)) - |T„(^)|2 = ±iC(^)i?„(E)e*"«"(^) ±isin(v^-7re„(E)) 

is real analytic with respect to a G M and hence Fj^\a) = C^^^ for all a G M. Also, 
it follows that the zeros of cos2(V^ — n^aiE)) — \Ta{E)\'^ (if any) are all of even order. 
Thus, we have that if E > belongs to one of the sets 5**^^^ then either 

cos^iVE-TtUE)) < \Ta{E)\\ (5.10) 
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or 

cos\VE - nUE)) > |T„(E)p, (5.11) 

for all a e M. Taking a = in dSTTt ) and using ^a=oiE) = and T„=o(^) = 1 gives that 
E = (vm)^ with some n eN. 

Now calculating the limit a ^ of Fj^\a) and taking into account that Ra{E) is not 
identically zero, we obtain that for E G 5^^^ 

sin =1= I sin Ve\ = 0, 

respectively. Thus, 

oo 

C |J[47rV,7r2(2n + l)% (5.12) 

n=0 

(5.13) 



C |J[7r2(2n + l)2,47r2(n + l)2]. 



n=0 



Now let us suppose that E belongs to the interior of the sets on the r.h.s. of ( |5.12| ). This 
implies that sin > if E G 5^+) and sin < if E G S'^'l Then using R^iE) = 
Ra{E)\ exp{i6a{E) — in^aiE)}, we calculate Fl^\a) (with the corresponding choice 
of sign) for small a: 



T- 



e-^<^^(E) + Oi\R^iE)f). 



2\sm{VE-nUE))\ 
Thus F^\a = 0) = and therefore F^\a) = for all a G R. This implies that 



-i sm{VE - 7r^a{E)) ± \J \R^{E)\^ - sm^{\fE - 7T^a{E)) = 

for all a G M. Hence Ra{E) = for all a G M, which contradicts the assumption E ^ S. 
For the sake of convenience we summarize some of the established properties of the 

sets 

Lemma 5.5 The sets S^"^^ are at most discrete. More precisely, 

C {(7^n)^ n G N}. 
IfE G then either f^T^j or ( pj7] ) holds. 



Now we define 

S= (5^+) U S^-^) n > : cos^(/E-7r^„(E)) < |T„(E)|2 for all a G M}, 
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which is at most discrete. 



Theorem 5.6 For almost all E > the upper Lyapunov exponent 7(-E') > almost 
surely. More precisely 7(-E') vanishes for E E S U S and almost surely nowhere else. 



Proof. We split the proof in several steps and start with the case E E S. Calculating 
the eigenvalues of Aa{E)Aa{E)* one easily finds that 



lUE)]' 



\TaiEW 



+ 



\Ta{EW 



\TaiEW 



1 > 1, 



where the norm is understood in the operator sense. Obviously, || Aq(_E) || = 1 iff Ra{E) 
0. Therefore, if ^ G 5 then 



log 



n a-.h(^) 



]=-n 



<log J] ||A„^,(^)(E)|| =0, 



]=-n 



and hence 7(-E) = 0. 

To proceed further we define the family of auxiliary periodic Hamiltonians 



3) 



(5.14) 



The spectrum of every H^°''> is purely absolute continuous and has a band structure. We 
recall (see e.g. [0) that the discriminant Aq,(£') of ( |5.14| ) is defined by Aa{E) = 
Mi(l) + M2(l), where ui{x) and U2{x) are solutions of H^^'^Ui = Eui with the initial 
data Mi(0) = 'U2(0) = 1, ^(0) = ^2(0) = 0. Aa{z) is an entire function of z E C. The 
real solutions of the inequality |Aq,(_E')| > 2 determine the gaps in the spectrum of H^°'\ 
whereas |Aa(-E)| < 2 implies that E belongs to the spectrum. 

We say that E is in a gap of for some a G M if there is 5 > such that {E — 
6, E + 5) n a{H^^'>) = 0. 

Keller proved that for E > 0, 



A„(E) 



2C0S(VE - TT^aiE)) 



Thus, E > is in a gap of iff cos'^{Ve - n^^E)) > |T„(E)p 



Remark: We sketch another proof of this fact based on the Ishii-Pastur-Kotani theorem 
]. If for some E > cos^{Ve — n^aiE)) > {TaiE)]"^, then one of the eigenvalues of 



A±(a) 



cos{VE - n^aiE)) cos^iVE - TV^aiE)) 



Km' 



- 1, 



(5.15) 
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is strictly larger than 1 in absolute value. Therefore, for the Lyapunov exponent corre- 
sponding to H^°''^ we have 

7„(E) = lim ^ log||A.(E)"+-+i 

m,n^oo n + m + i 

> lim ^ logmax|A±(a)r+™+^ 

m,n^oo 72 + m + 1 ± 

= log max I A-t (a) I > 0. 

Since the spectrum of H^^'^ is absolutely continuous, by the Ishii-Pastur-Kotani theorem E 
lies in a gap. Conversely, if cos^(V^ — -K^aiE)) < \Ta{E)\'^, then both eigenvalues ( |5.15D 
lie on the unit circle and the Lyapunov exponent vanishes, 

7a{E) = lim \ log ||A,(i?)"+'"+i = 0. 

Hence E belongs to the spectrum. 



We proceed with the proof of Theorem 5.6. Let us suppose that E & W \ S . Let 
Ge C SL(2; C) be the smallest closed subgroup which contains the support of ^ (recall 
that ip is the image of if under the map a Aa{E)). To establish that 'j{E) > for 
almost all £■ > almost surely we will also use the sufficient conditions of the Furstenberg 
theorem [[iT], Theorem A.II.4.1, Proposition A.II.4.3] (in the complex form), i.e. for almost 
all > the group G^; is not compact and for any x G CP^ the set {Ge ■ x} C CP^ has 
more than two elements. 

We start with the second condition in the Furstenberg theorem. First we calculate the 
(non normalized) eigenvectors of Aa{E), 



|T,(E)| I -zsiniVE - nUE)) ± \ cos^iVE - nUE)) - Km 



9t\E) 



In the case cos^{VE — TT^a{E)) = \Ta{E)\'^ both eigenvectors coincide and the corre- 
sponding generalized eigenvector (i.e. the solution of {Aa{E) — l)ga{E) = g^\E)) is 

(0,-if. 

By means of the bijectionp : CP^ — C U {oo}. 



P 



a;i ^ _ f ff, a;i ^ 0, 



X2 J [ OO, Xi = 0, 

we can identify the complex projective line CP^ and C U {oo}. Obviously, 

P{£\E)) = F^^\ 



[a] 



For fixed E the reflection amplitude Ra{E) as a function of « G supp(y9 has at most a 
discrete set of zeros. 

Actually (see [ pi] . Problem 6.4]) for any subgroup G C SL]r(2; C) one has the alterna- 
tives: 
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(i) G is finite, 

(ii) there is Q e GL(2; C) such that 

(iii) there is Q e GL(2; C) such that 

g-'Ggc|(^JJ ^, aeC\{0}, 6ec|, 

(iv) for any x e CP^ the set {G • x} c CP^ has more than two elements. 

Let us first suppose that G — Ge^s finite. Since Aq,(£^) is real analytic with respect to 
a this implies that Aa{E) is constant for all « G M and hence \Ta{E) \ does not depend on 
a. Therefore |T«(^)| = |T«=o(^)| = 1 and thus E e S. 

Consider the case (ii) for G = Ge- Let us suppose that for some Q the matrices 
Q~^Aa{E)Q are diagonal for all a in a subset of supp9? of positive measure. The existence 
of such Q implies that the eigenvectors of all Aa{E) are constant as elements of CP^ and 
thus 

p{g^^\E)) = C(±) = const (5.16) 

with C*^"*"^ 7^ C^~^ for all such a. But then these relations hold for all « G M by real 
analyticity. Thus E G ^^n^^"). Conversely, if E G S^n^W (such that cos^ = 1) 
and if in addition cos^(\/E' — 7r{a(£^)) ^ \Ta{E) f for all a in a subset of supp<^ of positive 
measure, then 



'^^^ pJ'iE)) Ptei-'W) )'^°L(2;C) 



does not depend on a in this set and Q~^ka{E)Q is diagonal for all these a. Hence for 
this choice of Q and by real analyticity A^Q is diagonal for all a G R. By the previous 
discussion (see Lemma 5.5) in case cos^(v^ — Ti^a{E)) < \Ta{E)\'^ with equality for a 
only on a set of measure zero the eigenvalues of Aq,^ (^) (E) lie on the unit circle. Therefore, 
since 

3 — 



we have || n,"!-™ A,^,(^)(S)|| = L Thus7(^) = 0. 

Consider the opposite case when cos^{\fE — n^aiE)) > |T„(£')p with equality for 
a only on a set of measure zero. In this case one of the eigenvalues X±{a) (say X+{a)) 
for almost all a G supp</? is larger than 1 in absolute value. Since E G 5'^+^ n S'^~'> the 
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eigenvalues of UT=-n^a,H{E) are given by UT=-nK{aj{uj)) and UT=-n^-{(^j{^)) 
and thus 

^ m 

l{E) > lim ^— _- ^ log|A+(a,M)| 

j=-n 

= E{log|A+(«,M)|}>0 

almost surely. Here we used the fact that {aj{uj)}j(zi is metrically transitive and the 
Birkhoff-Khintchin theorem. 

It remains to consider the case when E E S^^'^ fl 5^^^ and cos^(v^ — TT^aiE)) = 
|Ta(-E)P for almost all a in suppy^. By real analyticity this relation holds for all a; G M. 
But then A+(a) = A_(a) = ±1 and hence Q-^A^{E)Q = K{E) = ±1 for all a G M. 
Evaluating at « = gives Ta{E) = 1 for all a and thus E E S. We note that from the 
real analyticity with respect to a it follows that if there exists Q E GL(2; C) such that 
Q~^Aa{E)Q is diagonal for a on a set of positive measure, then Q^^ Aa{E)Q is diagonal 
for all aER. 

To conclude the discussion of case (ii) for G = Ge suppose now that there is Q G 
GL(2; C) such that for almost all a E suppv? 

0-A.(£)0 ''"J" ) (5.17) 

holds for some b{a) G C \ {0}. By real analyticity in a of Aa{E) such a relation holds for 
all q; G M. Taking traces gives 



C0S(V^ - 7T^a{E)) = 

for all a G M. Again by real analyticity this implies ^a{E) = const = ^a=oiE) = 0. In 
particular, ( |5.17D cannot hold if cos 7^ 0. In case cos ^/E = observe that for 6 7^ 



and 



with 



Q-'A^{E)Q^ 



i 
-i 



^ Ra{E) Ra{E) \ 

\ -t sin ^/E -i\To,{E)\ -t sin ^/E + t\Ta{E)\ J 



a 



Hence any matrix which diagonalizes Aa{E) is necessary of the form Qa ^ q ^ j ^^^^ 
ad ^ 0. This implies that 
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for almost all a G suppv?, where a{a) and d{a) are suitable nonvanishing functions. Writ- 
ing ( I5.18D explicitly gives 

qi = b{a){a{a) + d{a))Ra{E), 

q2 = —i{a{a)—d{a))Ra{E), 

gs = —ib{a){a{a) + d{a)) sin Ve — ib{a){a{a) — d{a))\Ta{E)\, 

94 = —{a{a) — d{a)) sin Ve — {a{a) + d{a))\Ta{E)\. 



From this it follows that for all a G supp^j 



Ra{E)q-i 
Ra{E)qA 



-iqi sin + q2\Ta{E)\, 



-iq2 sm 



^+-^AT^{E)\. 



(5.19) 
(5.20) 



From both ( |5.17[ ) and ( |5.20D it follows that b{a) ^ is real analytic with respect to « G M, 
and thus the relations ( |5.19[ ), ( |5.20[ ) hold for all a G M. Taking a = gives 

—iqi sin \fE + g2 = 0, 

77^ - iq2 sin/E = 0. 
6(0) 

The existence of a nontrivial solution in qi and q2 of this system implies that 6(0) = — 1. 
Thus, we obtain qi = —i sin \/Eq2 and qiq2 7^ 0. Inserting this in ( |5.19[ ) we obtain 
Ra{E)qs = g2(|T„(E)| - 1) for all a G R. This together with |T„(E)|2 + \Ra{E)\'^ = 1 
obviously implies that \Ra{E) \ is constant and thus Ra{E) = for all a G M. Thus E E S. 
This completes the discussion of the case (ii). 

Now consider the case (iii). As it is easy to see in this case that all Aa{E) have 
a common eigenvector, such that only one of the functions p{g^^\E)) is constant or 
p{g'i^\E)) = p{gt\E)) = const (and thus cos2(v^ - 7r^„(E)) = \T^{E)\^) for all 
a G suppv?. In the first case either E G S'(+) or ^ G S^'^ but E ^ S^+^ n S'-'\ In the 
second case the matrices Aa{E) are not diagonalizable and 

E G n |e > : cos2(v^-7rea(^)) = |T„(E)|2 for alia G m| 

= S^-^nf^E>0: cos2(/E-7r^„(E)) = |T^(£;)|2 for alia gm}. 

Conversely, if one of p{g^^ (E)) (say p{g^^ (E))) does not depend on a, then the matrix 



Q 



1 

,(+) 



Pig'cT'iE)) 



is such that 



Q-'Aa{E)Q 



A+(a) 






-1 



A„(a) 
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Let E e and cos^{y/E - iT^a{E)) < |T„(E)|2 for all a. Then all eigenvalues of 
all Aa{E) lie on the unit circle. We have 



A+K-(c.)) ^ 



n A„,H(i?)=g - n 

j=-n j=-n 

Obviously, there is 5 > such that 



(^) 



\Ro.{E)\ 



< 6 



\Tam 

for all a G suppv^. It is easy to see that if the numbers jSj are such that \/3j\ = 1, then 



n 



j=~n 



(3j bj_ 
(3, 



nr=.„/?. 



(5.21) 



with b satisfying the inequality |6| < J2f=-n I^j I- Thus, the norm of the matrix on the r.h.s. 
of ( |5.21| ) is less or equal a/2 + Therefore, 



n ^a.H(^) 



j=-n 



2x211/2 



< [2 + (n + m + 1)2^2] 



and thus 7(E) = 0. 

Let now E G 5*^+) and coi?{VE - 7r^Q(£')) > \Ta{E)\^ for almost a G suppv^. In 
this case either |A+(a;)| > 1 for almost all a G supp(/3 or |A+(q;)| < 1 for almost all 
a G suppv9. The eigenvalue of Y\^=~n^aj(u)){E) corresponding to g^\E) is given by 
UT=-nKiaji^))- Therefore 



log 



n ^".h(^) 



]=-n 



> ^log|A-,(a,M)| 

j=-n 



if |A+(q;)| > 1, or 



log 



>- 5^1og|A4a,H)| 

j=-n 



n ^-.h(^) 

j=-n 

if |A+(q;)| < L By the Birkhoff-Khintchin theorem we have 

^ m 

^— ^ log|A+(a,H)| =E{log|A+(a,H)|} 



lim - 

m,n^oo m -\- n -\- 



j=-n 



almost surely. If |A+(a;)| > 1 (< 1) for almost all a, we have E{log |A+(aj(t<;))|} > 
(< 0), and hence 



7(i?) = lim 



1 



m,n^oo m + n + 1 



log 



n ^-.h(^) 



j=-n 



> 
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almost surely. Similarly, in the case E G S^^"^ and cos'^{y/E — -K^aiE)) > \Ta{E)\'^ for 
almost all a E suppv? we have 'y{E) > almost surely. 

Thus, we have shown that the cases (i), (ii), (iii) occur iff G S U S^^^ U S^~\ 
Moreover, -f{E) = for E e S and 7(E) > almost surely for E G (5^+) U S^-'>) \ S. 
For those and only those E > 0, which do not belong to S' U 5''^+^ U S*^"^ the case (iv) 
occurs, and thus the second condition of the Furstenberg theorem is fulfilled. 

We turn to the first condition of the Furstenberg theorem. 

If is in a gap of H^°'^ for at least one a G supp(^, then one of the eigenvalues (|5.15p 



is strictly larger than 1 in absolute value. Therefore ||AQ,(i?)"|| > | max± A-|-(a)|" ^ 00 as 
n — > cx), and hence G^; is not compact. 

Now suppose that E is not in a gap of H^'^'^ for all a G suppy^, i.e. cos^{VE — 
7r^a{E)) < |Tq,(£')P for such a. In this case the eigenvalues \±{a) lie on the unit circle. 

Since suppy^ contains at least two a's such that the corresponding Aq,'s have no com- 
mon eigenvectors, by the assumed continuity of (p we can select them in such a way that 
either cos2(/E-7r^„(E)) < \T^{E)\^ or cos\y/E - n^^{E)) = |r„(E)|2 for both a's. In 
the first case we can apply the arguments of pB] , to construct a matrix, which belongs 
to G^; and has an eigenvalue strictly larger than 1. To treat the second case we consider 
two matrices 



b* a* 



with no common eigenvectors, satisfying |Reaj| = 1 for both j = 1,2. Since Af, G 
SL(2; C) one has |Imaj| = \bj\. The matrices Mj are not diagonalizable, because to the 
eigenvalue \j = Reuj corresponds the only eigenvector gj = (ibj,liaaj)'^ . Further we 
consider the matrices Mi, M2, which are defined as follows: Mj = M|, j = 1,2 if 

ReaiImaiRea2lma2 < 0; Mi = Mf, M2 = M2"^ if ReaiImaiRea2lma2 > 0. These 
matrices have the form 



M, 



aj bj 
b* a; 



such that Reuj = (Rea-,)^ = 1 and Imailma2 < 0. We can write the matrices Mj in the 
following form 



Mi = I + zIma,A^. 



with 



1 e*'^^' 
= ( -e-'^J -1 

where e^^^ = —ibj/lmaj = —ibj/lmaj. Obviously, Nj = 0. Since by assumption the 
eigenvectors of the matrices Mj are linearly independent, we have that ^ ^2- Now we 
calculate 

, ^ ^ 

tr(MiM2) = 2 - 41mailma2 sin^ — > 2. 



5 Lyapunov Exponent 



45 



Therefore, one of the eigenvalues of M1M2 is strictly larger than 1. 

Thus, by the Furstenberg theorem it follows that 'j{E) > almost surely for E ^ 
S U S^~^^ U S^~\ This completes the proof of the theorem. ■ 

Let (■, ■) denote the inner product in C^. We note that 

^(-'^,™)(^. ^) = ±1 arg (e±, A(-"''")(E; uj)e±) 
n 

with e+ = (1, 0)^, e_ = (0, 1)^. From Q it follows that 

/ m 



]=-n 
m 



\ j=-n 

Obviously, UEe± = e'^^^-^e±. Therefore the r.h.s. of ( |5.22D equals 
Thus we obtain 



j=-n 



TT TT n,m^oo n + TTl 



\ j=~n / 



Since ^{E) = No{E) - N{E) = VE/tt- N{E), it follows that 



^-^argle±, J| Aa,iu.){E)eA 

\ j=-n / 



(5.22) 



N{E) = T- lim - 

TT n,m-*oo n -\- m 

This representation is similar to the definition of the density of states through the rotation 
number of fundamental solution of the Schrodinger equation . 

The representations (|577|), ( |5TP| ) can also be rewritten in a similar form, 

7(i?) = lim — ^^log|(e±,A(-"''")(^))| 



m,n^oo m + n + 
1 



lim - 

m.n^oo m + n + 1 



log 



]=-n 



Example 1. Here we consider the Hamiltonian ( PTTj ) where / is (formally) replaced by 
the Dirac 5-function. In this case the transmission and reflection amplitudes are given by 



T„(E) 

Ra{E) 



1 + 



2/E 



a 



—i- 



1 + 



la 



i4E \ 2VE 
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Therefore 5 = and 



2 \/ E 

F^f ^ (a) = ( sin + cos 



2\/E // ^ a r-\ 
±i-^W cosV^+ — ^sinv^E - 1, 



cos^(v^-7rea(^)) / /^^ « • /^^ 



2 



Thus 5 = {(vr/c)^, k E N}. From Theorem 5.6 it follows that 'j{E) is positive on (0, oo) 
except for the set S. That 7(-E) = iff E = = (vrA;)^, G Z was proved by Ishii 
(see also \^). Kirsch and Nitzschner ^] proved that A^(Efc) = A^o(^fc)- Here we 
reconsider these facts once more. 

The matrices Aj{E) can be calculated explicitly: 



A,{E)=I + ^j^A,{E), 



where 



ME) 



For E = Ek,keN 



Aj{Ek) = A = 
Obviously A is nilpotent, i.e. = 0. Therefore 



I ^-2iVEj 



1 1 

-1 -1 



K^-'-\E,-u)= n A.™= n v+Ywv 

i' — i' — _nV w k / 



From this it follows that 



j=-n 3=-n 

m 



2 V il^i. 



2 



(e±,A(-"'™)(i?,;a;)e±) = 1 ± ^ 



2v 

V 1^ j=-n 



and hence 

arg (e±, A^""''"^(Efc;u;)e±) = iarctan 



e±,A(-"'"^)(E,;c.)e±)| = l+_ 



\j=-n / 
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Since {aj{Lo)}j^z is metrically transitive, by the Birkhoff - Khintchin ergodic theorem we 
have 

lim } aj(uj) = lim > aoiTjCu) = Ejaj. 

n,m^oo m + n + 1 ^-^ n,m-»oo m + 77, + 1 

j=-n j=-n 

Therefore, the sum J2f=-n'^ji^) limit m,n ^ oo increases not faster than (n + 

m + 1). Hence 

C{Ek) = - lim / arg(e±,A(-"'-)(E,;a;)e±) =0, 

7(Efc)= lim / log|(e±,A(-"''")(E,;cu)e±)| =0. 

6 Analytic Continuation and the Thouless Formula 

In this Section we show that the analyticity of w(E) = —'j{E) + i7rN(E) in the upper 
complex -plane C+ and the Thouless formula ( p3| ) are a direct consequence of the fact 
that the functions logT(-"''")(E) = log |T(-"''")(E)| + i7r^(-"''")(E) are analytic in C+ 
for every n, m G N. 

By Lemma 2.3 we have 



logT(-"'™)(2) 



E-z 



Lemma 6.1 Let Eq = max-t {|q;±|} /o with /o = sup Then there is a constant 

C independent ofE and n,m E N such that 

|ei-"'™)(i5)|<-^(n + m + l) 



for all E > Eq. 

Proof. By the monotonicity of the spectral shift function we have 

|^(-n,™)(i^)| < ^^E- H, + E,x^--'^\ Ho), 

where x^""'*") is the characteristic function of the interval [— 72 — 1/2, m + 1/2]. Calculating 
({E; Ho + EoX^""'™\ ^o) explicitly we find 



aE; Ho + i?oX^-"'"^\ Ho) = ^{n + m + l) 

TT 



N 1, f 2E-Eo , , , ,A 

[n + m + 1) Arctan — — — tan(-\/ — Eo[n + m + 1)) , 

TT TT \1\/E-JE - Eq J 
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where Arctan is the mukivalued arctan function such that Arctan(C tan(a;)) is continuous 
and nondecreasing with respect to x. Since 

2E — En 

" > 1 



24E^E - Eq 

for all E > EQ,ii follows that 

f 2E-Eo \ , 

Arctan — — tan x > Arctanltanx) = x. 

\2^/E^E-Eo J 



Therefore 



m Ho + EoX^-'-\ Ho) < v^-v^^^ (, + ^ + 1^ 

TT 



E - -JE-Eo < ^° 



Obviously, 



'E 

for all E > Eo, thus proving the lemma. ■ 



lim sup {n + m + 1)-' I '^'^ , '^ dE = 0. 



Now to study the limit ra, m ^ oo we can use the theorem on the continuity of the 
Stieltjes transform (see e.g. Appendix A]). The applicability of this theorem is guar- 
anteed by Lemma 6.1, from which it follows that 

Jr l + \E\ 

and 

hm fiur) (n. -\- m. -\- 1)^^ I 

1^1 

Therefore, since {n + m + -> ^{E) for all E e M and P-almost all G fi, 

and ^{E) is nonrandom, we obtain that {n + m + 1)^^ log Ti""'™-* (i?) converges for all 
z E C with Imz > and P-almost allcu E to some deterministic limit W{z), which is 
given by 

Jr — z 

and therefore is analytic for Imz > 0. 

Since ^ (E) is continuous (moreover Holder continuous) by the Sokhotski-Plemelj for- 
mula we have that W{E + iO) exists for all G M and 



lmW{E + iO) =mi{E). 
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As in the proof of Lemma 2.3, since ^{E) is continuous and is of bounded variation, we 
get 



ReW{E + iO)= I \og\\ - E\d^{\). 
On the other hand, by Lemma 2.3 we have 

Jr 

By Theorem 5. 1 we have that for every fixed E > 



(6.2) 



(6.3) 



-liE) 



lim / loK |A — E\ 



7i,m^oo /,n n + m + 1 



(6.4) 



almost surely. Now we prove that for Lebesgue almost all ii^ G M there are subsequences 
Uj, rrij tending to infinity and such that for P-almost all cu 



lim / log |A — E\ 

j^oo Jjg rij +mj + 1 



\og\X-EmX). 



(6.5) 



Thus, from (|]2|) and (^ it will follow that ReW{E + iO) = --f{E) and therefore 

- / log|A-i?|d^(A) 



for almost all E > 0. 

The arguments used below are very similar to those of Pastur and Figotin pSj Theorem 
11.6]. Consider the functions 



T, 



(—n,m) 



t{E) = [ \og\X-EmX)= [ log 
Jr Jr 

d^(.~^'"^)(X) 



X-E 



log X-E 



n + 171 + 1 

Let us fix some interval K cM. and consider 

tB{\) 





X-E 


/ log 


IB 


i-E 



i-E 
log 

dE, 



d^X), 



X-E 



i-E 



rf^(-'^'™)(A) 
n + m + 1 



where B E B{K) (the set of all Borel subsets in K). The family of functions {tsiX), B G 
B{K)} is uniformly bounded and equicontinuous on any bounded interval and 

sup |tB(A)| <C(l + |A|)-i 

B&B{K) 



for some C > 0. Also we have 



T{E)dE = / tB{X)d^{X). 
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Since the family tsiX) is uniformly bounded and equicontinuous and since (n + m + 
1)~M.^(~'^''")(A; iu) converges vaguely to d^{X) almost surely, it follows that 



lim sup 



{tI-'-\E) - t{E)) dE 



0. (6.6) 



On the other hand one has 



/ \rl~^^'^\E) - T{E)\dE 

JK 

[ {ri~-'^-\E) - r{E)) dE 

-/ {Ti--'-\E) - r{E)) dE 

[ [T^--'^\E)~T{E))dE . 

J B 



< 2 sup 

BeB{K) 



From ( p^ it follows that 

lim / |r^-"''")(E) -r(E)|dE = 0. 

n,m— >oo J 



Denote by Vli the subset of Vl such that P(r2i) = 1 and for which (n + m + 1)~^ 
d$,^~'^'"^\X; Lu) converges to d$,{X). Taking an arbitrary G fii we can choose subse- 
quences rij, nij tending to infinity and such that 

lim |r^-"^'"^)(E) - t{E)\ = (6.7) 



for Lebesgue almost all E E K. Since {n + m + l)~^c/^(~"'™)(A; to) — > d^(\) for almost 
all a; e fi, relation ( p^ remains valid for almost every u) E il. Thus (^3[) is proven. 
We have shown that 

lmW{E + iO) = mi{E) 
ReW{E + zO) = -7(E) 

for almost all > 0. Therefore, W{z) is an analytic continuation of —'y{E) + m^{E) 
from the real semiaxis. 

Let us consider the function w{z) = W{z) + ^0(2:), where 

Wo{z) = -7o(^) + iTrNo{z) = 

with 7o(-E') being the Lyapunov exponent for Hq = —d^/dx^. It is obviously analytic in 
the upper half-plane \m.E > and 

w{E + iO) = -i{E) + i7TN{E), E > 0. 

From ( |0| ) and from the fact that ^ (E) is given by the difference of two nonnegative 
functions No{E) and N{E), it follows that W{z) is the difference of two Nevanlinna func- 
tions. Since wo{z) is a Nevanlinna function, so is w{z). Also from ( ]6.5D it follows that the 
Thouless formula (Ob holds. 
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7 Some Extensions 

Let V{x) be some real-valued uniformly locally integrable function. For simplicity we 
suppose that V{x) is uniformly bounded, i.e. < Vo,Vo > 0. The last assumption 

can be weakend but we do not go into details here. Let {yjjjez be a sequence of real 
numbers such that yj — > ±oo as j — > ±oo and yj < yj+i for all j e Z. We suppose that 
{%}iez is such that all the differences yj+i — yj are finite (but not necessarily uniformly 
bounded). Let Xji^) be the characteristic function of the interval [yj,yj+i]. We denote 
Vj{x) = V{x)xj{x) such that 

m 
j=-n 

tends to V{x) as m, n ^ oo. By the above assumption one has |V^(a;)| < VoXjix). 

Let H, Hj and denote the Hamiltonians with domains of definition being the 

Sobolev space (2,2) (^^^ 

H ^Ho + V, Hj ^Ho + Vj, ^Ho + V^^"'^'"^). 

Let ^j{E), Tj{E) and be the spectral shift function and transmis- 

sion amplitude for the pairs (Hj, Ho) and Hq), respectively. Also as above we 

denote 

gr"'™)(E) = e^-"''")(E)±i. 

By Corollary 3.4 and Theorem 5.2 we have 

for every k E Z such that — n < A; < m. By the monotonicity and superadditivity proper- 
ties of the spectral shift function (Corollary 3.4) we also have 

{ym - y-n)-'l^r''^\E) > {ym ' IZ-n)"' [m ^0 " , ^^o) + l] 

> inf (y„ - [i{E- Ho - VoX[y.„,y^], Ho) + l] 

m,n 

= lim (y^ - [^(E; Ho - VoX[y_„,yra], Ho) + l] 

m,n—^oo 

= -[max{0,E + Vo)Y/yrr 

for all E eR. Similarly, 

(yrn - y.n)~'f-''''^\E) < {ym - y-nV [i{E- Ho + VoX[y.^,y^].Ho) - 1] 

< [max(0,^- K))]^^V^- 
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Also we have that |T(-"''")(E)| < 1. Therefore, if Hq + l/C""''") and Hq + yC-^+^.^+fc) 
are unitary equivalent from the known property of subadditive functions (see e.g. [ p3 , 
Theorem 6.6.1]) the existence of the limits 



m 

1t{E) 



lim , 

log|T(o.™)(^)| 



— lim 

m— »oo 



log|T(-"'0)(E)| 
lim ■ ■ , 

\y~n\ 

log|T(-"''^)(E)| 
lim 

m,n^oo - y_„ 



follows. Clearly 7^, 7^,, and 7^ may be unequal. Theorems 4.3 and 5.3 apply also to this 
case. Thus we have again that ^(E) = No{E) - N{E) for all E e M and 7^(E) = -f^{E), 
7t(-E) = liE) for all positive E, where 



iHe) 

l{E) 



lim — log||0(O,a;;E)||, 
lim ^\og\\(t){-x,x]E)\\ 



are the upper Lyapunov exponents. Here x'; i?) denotes the fundamental matrix of the 
Schrodinger equation with the potential V . Further we can again prove that 



l{E) 
N{E) 



lim log 

m,n^oo y^ - y_„ 



n ^^■(^) 



]=~n 



1 



=p— lim - 



^ arg I e±, JJ Aj(^)e± j 

V j=-n J 



where 



K,{E) 



M£)piv^{%+i-2%+W-i)/2 



.M^p-«v^fe+i-2%+%~i)/2 



TjiE)' 



»VE(yj+l-!/j_l)/2 



and Tj{E), Rj{E), Lj{E) are transmission and reflection amplitudes corresponding to the 
potential Vj{- — Uj). 



Appendix A 

Here we prove that the expression in the braces on the r.h.s. of ( |5.3D is bounded from below 
by the constant C{E) = AE/{1 + E"^). For brevity we set B = J I - |Ti'"|2, a = smOi, 
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and b — sin 62. Then the expression at hand can be written as follows 

F{a,b,E,B) = 2^2(1 -a^) + 2(1 -62) 

+hBa-bf + E{Ba + bf. 
E 

We show that F{a, b) > C{E) for all > independently of a, 6 e [-1, 1] and S e [0, 1]. 
First we note that 

F{a,b,E,B) = F{-a,-b,E,B) 
= F{a,-b,l/E,B) = F{-a,b,l/E,B). 

Since C{E) — C{l/E) it therefore suffices to consider the case a, 6 e [0, 1]. Since 

d 

777-F(a, b, E, B) = 25(2 - 2a^ + a^/E + a^E) + 2ab{E - 1/E) 
oB 

is nonnegative whenever E > 1 we have 

F(a, 6, E, B) > F(a, b; E, 0) = 2(1 - ft^) + b^/E + Eb^ > 2 > C{E) 

for such E. Here and in what follows we use the estimate E + 1/E > 2 for all > 0. 
Now take the case < < 1. Then for fixed E, a and b the function 

^F{a,b,E,B) 

has exactly one zero as a function of B in the interval [0, 00) at 

-ab{E-l/E) 
° ~ (E + l/E-2)a2 + 2' 

which is a minimum for F{a, b, E, •). Hence 



F{a, b, E, B) > F{a, b, E, Bq) 

a^b^{E-l/E)^ 
~ {E + 1/E -2)a? + 2 
=: G{a,b,E). 



2 + {l/E + E -2)b^ 



It is easy to see that 



d ^, , Aab'^iE - 11 Ef 



Therefore 



G{a,b,E) > G{l,b,E)^2 + 2b 
> G{1,1,E) = 2 + 2 



22- 1/E-E 



E+l/E 
2-1/E-E AE 



E + l/E 1 + E^' 
thus proving the claim. 
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